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We study the question whether matrix models obtained in the zero volume 
limit of fd Yang-Mills theories can describe large N QCD strings. The matrix 
model we use is a variant of the Eguchi-Kawai model in terms of Hermitian 
matrices, but without any twists or quenching. This model was originally 
proposed as a toy model of the IIB matrix model. In contrast to common 
expectations, we do observe the area law for Wilson loops in a significant 
range of scale of the loop area. Numerical simulations show that this range 
is stable as N increases up to 768, which strongly suggests that it persists in 
the large N limit. Hence the equivalence to QCD strings may hold for length 
scales inside a finite regime. 
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1 Eguchi-Kawai equivalence 

In quantum field theories and systems in statistical mechanics, the number of 
internal degrees of freedom often enters as a free parameter, which we denote 
by N. In many cases, the dynamics of the system simplifies significantly as 
N becomes large. A general discussion is given for instance in Ref. (ij. 

This effect motivated in particular the model suggested by Eguchi and 
Kawai a long time ago ||2|. Their point of departure was standard U(N) 
lattice gauge theory. Based on the factorization of correlation functions at 
N —> oo, they suggested that the model should be equivalent to its dimen¬ 
sional reduction to one point. Then all the link variables are replaced by 
U Xitl —> I/p, and the plaquette action reduces to 

Sek = ~N(5 J2 Tt(U,U u UIuD . (1) 

As their argument for the Eguchi-Kawai equivalence to ordinary gauge the¬ 
ory, the authors showed that the Schwinger-Dyson equations remain unal¬ 
tered, which also implies the invariance of the Wilson loops. However, their 
derivation implicitly assumed that the [U(l)] d symmetry of the phases is not 
spontaneously broken. This is correct at strong coupling in dimension > 2 
and for all couplings in two dimensions. On the other hand, it turned out 
that at weak coupling in dimension > 2, spontaneous symmetry breaking 
does set in. Numerical simulations revealed that the N eigenvalues of U M 
tend to be almost equal [|J . Hence this attempt to prove general equivalence 
to standard lattice gauge theory failed, and therefore the equivalence to the 
continuum gauge theory is uncertain as well. 

In order to avoid this problem, soon afterwards modified versions of the 
Eguchi-Kawai model were introduced, in particular the “quenched Eguchi- 
Kawai model” U and the “twisted Eguchi-Kawai model” ||, where the phase 
symmetries are not spontaneously broken. Ref. || also contains numerical 
results on the quenched Eguchi-Kawai model, and Refs. &1 present a per¬ 
turbative consideration. For numerical studies and generalizations of the 
twisted Eguchi-Kawai model, see for instance Ref. 0. 

More recently, Ishibashi, Kawai, Kitazawa and Tsuchiya worked out a 
large N reduced matrix model of lOd super Yang-Mills theory [j|] (for a re¬ 
view, see Ref. 0). The authors presented some evidence for its relation to 
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type IIB superstring theory and conjectured that it could provide a con¬ 
structive (non-perturbative) definition of string theory. In particular, they 
showed that the matrix model action reduces to the Green-Schwarz-Schild 
action in a certain semiclassical limit []. Moreover, their IIB matrix model 
is shown to describe D-branes and their interactions correctly, which was 
further elaborated in Refs. [11 . 

Supersymmetry seems to protect the [I/(l)] d symmetry [§], but it turns 
out to be marginally broken (l^|. It is completely broken in the simplified 


bosonic IIB matrix model, which is obtained by omitting the fermions. Nev¬ 
ertheless, the bosonic model has been recognized as a well-defined model 
13], and it has attracted considerable interest in the literature. Both, the 


bosonic and the supersymmetric model can be formally obtained from the 
zero-volume limit of large N gauge theories in the continuum, and hence they 
can be viewed as variants of the original Eguchi-Kawai model. The finite¬ 
ness of the partition function becomes a non-trivial issue, since the matrices 
to be integrated are Hermitian. The condition for the finiteness (N > 2, 
d = 4, 6,10 for the supersymmetric case, and N > for the bosonic 

case) conjectured in Ref. ||T3|, [L4[] has been proved recently |T5|. However, the 


validity of Eguchi-Kawai equivalence is still an open question. 

In particular, the 4d case of this Hermitian matrix model is accessible for 
numerical studies. This is true even for its supersymmetric form, because the 
fermionic determinant is real positive in d — 4 |L6|, as suspected earlier [|T|] . [] 
Large N factorization is known to hold in the 4d IIB matrix model, both, for 
the bosonic jTT| and for the supersymmetric case |TB[. As a stronger criterion 
for Eguchi-Kawai equivalence to ordinary gauge theory, one may investigate 
the validity of the area law for the Wilson one-point function. Simulations 
for N up to 48 show a finite range of the scale, where the area law seems 
to hold 0. [] For the SUSY case, larger N can hardly be simulated on 
present computers, but in the bosonic case it is possible to go far beyond. [[ 
The goal of this work is to investigate the fate of the finite area law window 


4 A more in depth study of this correspondence has been reported in Refs. |To| . 

5 The bosonic model has also been simulated in higher dimensions. In particular, Ref. 
0 presents results up to d = 20, N = 32, which are compared to a 1/d expansion, and 
Ref. |I^| arrives at d = 10, N = 128. Also the SUSY model without fermionic phase factor 
was simulated in d = 6 and d = 10 |l9| . 

6 A similar behavior was also observed in the lOd bosonic model ||l8f . 

7 The computational effort grows like N 4 5 6 7 in the SUSY case, but only like N 3 in the 
bosonic model ||l^|. 
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as N becomes really large. We increase N up to 768, and observe that the 
scaling window does indeed stabilize; this provides strong evidence that it 
does neither shrink to zero nor extend to a larger regime in the large N limit. 
We also confirm for N up to 768 that the correlation functions relevant for 
QCD and string theory remain finite in the large N limit if one keeps the 
product g 2 N fixed [[Hi . 

In Section 2 we describe the model we are studying, and the properties 
which are relevant in this context. In Section 3 we recall the notion of Wilson 
loops in the framework of matrix models and their scaling behavior. Then 
we present our results for the area law at N — 48, 128, 256, 512 and 768 for 
the bosonic 4d IIB matrix model. In Section 4 we discuss the conclusions 
from these results, which were mentioned before in Refs. [EO]. 


2 The model 


In 1997, Ishibashi, Kawai, Kitazawa and Tsuchiya suggested a model of su¬ 
persymmetric matrices, which is a candidate for a non-perturbative definition 
of IIB superstring theory ||. Here we are interested in the bosonic variant 
of this model, which is obtained if one simply drops the fermions by hand. 
The resulting bosonic IIB matrix model is given by the partition function 


Z = J dA exp (—S'[.A]) , (2) 

SM = -4. Tr([^,^] 2 ), 

where A M are traceless Hermitian N x N matrices, and in the d dimensional 
version of this model, /j runs from 1 to d. The original model was formulated 
in d = 10, but here we consider d — 4. The constraint of tracelessness is 
incorporated in the measure, 


dA= n 


M=1 


N N 

IK'®* (A^ij dIm(H At )jj} n w d ew 

i>j i=l i= 1 


( 3 ) 


Note that this model is well-defined without any cutoff for N > [R 

Hence its only parameter g is a pure scaling parameter, rather than a coupling 
constant. It can simply be absorbed by rescaling the variables as 


\fg X/i i 


( 4 ) 
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where the matrices X^ are dimensionless. 

This model is invariant under Euclidean rotations of the A^, and it has 
the SU(N) symmetry 

A , u VA^ , V e SU(N) . (5) 


This 4d bosonic IIB matrix model is closely related to the original Eguchi- 
Kawai model in the weak coupling limit (for details, see e.g. Ref. [17]). Due 
to the [£/(l)] d spontaneous symmetry breaking the eigenvalues of the model 
(jl]) collapse to a point and we obtain 


Sek ~ —yNfl E Tr {[A,, A u } 2 ) + 0(A 6 ), 

^v=l 


( 6 ) 


where U fl = e lAl \ Thus we can identify f3 = l/(2g 2 N ) and the two models 
agree in the limit /3 —> oo. 

Moreover, the large N limit of the bosonic model is expected to be equi¬ 
valent to the ’t Hooft limit of the corresponding Yang-Mills theory with 
a certain ’t Hooft coupling constant Aym and a momentum cutoff A (17 


In this correspondence, however, the ’t Hooft coupling constant is related 
to the momentum cutoff canonically, and therefore one cannot tune Aym 
appropriately as one sends A —> oo. 

In the next Section we are going to take a new look at the possibility of 
the Eguchi-Kawai equivalence of the bosonic model to ordinary gauge theory 
at large N. After all, the presence of spontaneous symmetry breaking does 
not rule out that Eguchi-Kawai equivalence could be realized for some range 
of length scales. As a powerful criterion to test this possibility, we study how 
far the Wilson loops follow an area law. 


3 The area law 

In analogy to ordinary gauge theory, Polyakov lines and Wilson loops are 
introduced as 


m = 

—Tr[exp(7fcAA)] , and 

(7) 

W(k ) = 

—Tr[exp(7fcAA) exp (ikX 2 ) exp(-ikXi) exp(— ikX 2 )} , 

(8) 
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respectively. For convenience we insert particular components of the dimen¬ 
sionless matrices X fJ , but their choice is irrelevant. In the interpretation of 
large N reduced matrix models as string theory, Wilson loops play the role 
of string creation or annihilation operators |21| . The parameter k represents 
the momentum distribution carried by each section of the string. Its physical 
(dimensionful) counterpart is given by fc p h ys = k/^/g. Correlation functions 
of these operators have a direct physical meaning, and hence the existence of 
a non-trivial large -N limit is absolutely crucial. It has been observed to exist 
if one tunes g oc 1 /\/iV, in the bosonic and in the supersymmetric case Jl6 . 
In what follows, we set g = \J48>/N following the convention in Ref. ||16|| . In 
the Eguchi-Kawai equivalence, (W(k)) corresponds to the Wilson loop in the 
large N gauge theory, and the parameter fc p h ys represents the linear extent 
of the Wilson loop in a physical scale. 

In Fig. |l| we illustrate the large N scaling further by showing the Polyakov 
line (P) as a function of k p h ys , with N = 48, 128, 256, 512 and 768. Scaling is 
confirmed in an impressive way. We also observe that the fc p h ys -dependence of 
the scaling function is in qualitative agreement with an analytical prediction 
by Oda and Sugino, which is based on a Gaussian approximation (to the 
second order) in the large N limit |22|. The observation that ( P) does not 


vanish at finite values of fc p h ys confirms that the [[/(l)] 4 symmetry is broken, 
as we mentioned in Section 1. 

Fig. || is the corresponding plot for the Wilson loop (IF), also as a function 
of k phys . We see that large N scaling is confirmed again, hence we can 
take a rather safe extrapolation to N —> oo. Again we show the Gaussian 
approximation (to the Erst order) |22|] for comparison. 

In Fig. [| we show the Wilson loop logarithmically against k phys , and 
we see a window, ranging approximately from k phys & 0.5... 1, which is in 
agreement with the area law. This area law regime turns out to be manifestly 
stable as N becomes large. Therefore, we do observe a finite range of scale 
where Eguchi-Kawai equivalence is in business. 

In Fig. f| we extend this plot up to quite large momenta in order to check 
if perhaps a further area law regime shows up. We observe an oscillating 
behavior, which also persists at large N. These oscillations might include 
additional, small area law regimes, but the relevant interval with this respect 
is clearly the Erst one, shown in Fig. [|. We also looked for a regime where 
the perimeter law holds, but we observed that there is no such interval with 
a size comparable to the leading area law regime. 
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Figure 1: The Polyakov line ( P) plotted against the dimensionful momentum 
k p hys = k/y/g. Numerical results for N = 48... 768 confirm the large N 
scaling, as well as qualitative agreement with the Gaussian prediction. 

4 Conclusions 

In view of the literature on the Eguchi-Kawai model, it comes as a surprise 
that the area law does hold in a significant range of scale for the bosonic IIB 
matrix model. It implies that Eguchi-Kawai equivalence may hold in this 
regime, even in the absence of twists or other mechanisms to preserve the 
phase symmetry. 

This observation opens up interesting perspectives, with respect to both, 
the theoretical background and applications. It is now motivated to recon¬ 
sider the possibility of Eguchi-Kawai equivalence also in the original Eguchi- 
Kawai model, or further variants thereof, which do not involve a symmetry 
preserving mechanism. In the matrix model studied here one could still in¬ 
vestigate whether the area law regime extends to infinity if we modify the 
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Figure 2: The Wilson loop (W) plotted against the dimensionful momentum 
^phys = k/y/g. Numerical results for N = 48... 768 confirm the large N 
scaling. Again the Gaussian approximation is shown for comparison. 

model by the quenching procedure |6j. This question is non-trivial given 
that the Hermitian matrix model has no direct link to lattice gauge theory, 
and its Eguchi-Kawai equivalence has been discussed only at the perturbative 
level [fj]. Such a study may also shed light on the string-theoretical meaning 
of the quenching procedure. 

We hope that our results presented here provide new insight towards a 
string description of large N QCD. 
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Figure 3: The Wilson loop plotted logarithmically against k 2 hys = k 2 / g. The 
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